I. Key to understand direct and indirect effects of climate change
  A. Particularly important wrt foundation species
II. Kelp forest paradigm is envt -> kelp -> community
  A. Some direct envt effect, but variable
III. New England kelp morphology very different, different ecology
  A. Levin et al indications that kelp is not only foundation species here
  B. Look at Ecklonia systems
  C. Look, we’re New England. Whatever, California.
  D. Transition with talking about sugar kelp as dominant
IV. Sugar kelp is sensitive to temperature via multiple mechanisms
  A. Temp up, kelp should go down from expts
  B. Review everything from Gerard to killer Nova Scotia work
V. But does it blend in the field?
  A. To answer if temp drives kelp and community, we put together large-scale survey
  B. KEEN intro, and how it will answer questions.
  C. Leave the reader wanting more by saying we show that our system is totes different
We collected temperature from the nearest offshore buoy to each site. While these temperatures were likely not the conditions experienced precisely at the site, they should at least correlate to both regional and temporal trends experienced by kelps. As kelp reproduction and initial growth occurs in the winter and early spring, we chose to examine minimum winter temperature as a reflection of conditions that could influence both kelp abundance as well as influence the community as a whole. note to self - this justification either needs beefing up, referencing - ha, or I need to systematically test a few other summarizations of temperature. They’re all pretty correlated anyway –Note to co-authors, I also attempted this briefly with nearest NOAA monitoring station, as this pulled stations within the Narraganset and Boston Harbor, but those stations differed by upwards of 8C on minima from others, and not in directions you would suspect, leading me to think that offshore is the more appropriate way to go
While some sites had benthic loggers deployed during the period of observation, buoy data enabled a more comprehensive analysis for two reasons. First, this meant that the first year sampled had to be thrown out due to no preceeding winter temperatures recorded. This filtering was particularly problematic as a few sites only have one year of data at present. Second, For analyses to determine a proper causal linkage (see data analysis below), we required temperature data during the duration of sampling for the entire data set - e.g., even for a site sampled in 2016, our models required data from 2013-2017.
To examine the relationship between temperature and Saccharina latissima abundance, we first took a naive look at the relationship between winter minimum temperature at kelp abundance by fitting a generalized linear mixed model with a lognormal error and log link. Percent rocky substrate was included as a covariate, and we included a random effect of site. To evaluate the model, we include a chi square log ratio test of the model as well as tests of the coefficients.
| Chisq | Df | Pr(>Chisq) | |
|---|---|---|---|
| WINTER_MIN_SEA_SURFACE_TEMPERATURE | 7.497 | 1 | 0.006 |
| PERCENT_ROCK | 1.055 | 1 | 0.304 |
| effect | group | term | estimate | std.error | statistic | p.value |
|---|---|---|---|---|---|---|
| fixed | fixed | (Intercept) | 1.644 | 0.708 | 2.323 | 0.020 |
| fixed | fixed | WINTER_MIN_SEA_SURFACE_TEMPERATURE | -0.146 | 0.053 | -2.738 | 0.006 |
| fixed | fixed | PERCENT_ROCK | 0.004 | 0.004 | 1.027 | 0.304 |
| ran_pars | SITE | sd_(Intercept) | 1.914 | NA | NA | NA |
| ran_pars | Residual | sd_Observation | 3.851 | NA | NA | NA |
While this model seems to demonstrate a negative relationship between temperature and sugar kelp abundance, given spatial gradients of temperature that correlate with a wide variety of other drivers, it’s difficult to say whether this represents a true relationship. To assess the importance of temperature with greater confidence, we used a hierarchical model with group mean centering [@REF TBA]. By calculating a centered mean temperature from 2014-2016, an annual deviation from site-level mean, and including both in the model along with a random effect of site, we’re able to 1) derive a relationship driven by between-year variation in temperature without the confusion of spatial correlates and 2) evaluate the relationship due to spatial differences in thermal regime via use of a hierarchical predictor alongside a site-level random effect. In essence, we partition the spatial and temporal component of temperature’s effects, and the deviation from site average coefficient addresses the impact of temporal change in temperature at a site.
| Chisq | Df | Pr(>Chisq) | |
|---|---|---|---|
| group_temp | 29.795 | 1 | 0.000 |
| temp_group_cent | 6.083 | 1 | 0.014 |
| PERCENT_ROCK | 1.299 | 1 | 0.254 |
| effect | group | term | estimate | std.error | statistic | p.value |
|---|---|---|---|---|---|---|
| fixed | fixed | (Intercept) | 3.230 | 0.336 | 9.619 | 0.000 |
| fixed | fixed | group_temp | -0.699 | 0.128 | -5.458 | 0.000 |
| fixed | fixed | temp_group_cent | -0.118 | 0.048 | -2.466 | 0.014 |
| fixed | fixed | PERCENT_ROCK | 0.004 | 0.003 | 1.140 | 0.254 |
| ran_pars | SITE | sd_(Intercept) | 0.000 | NA | NA | NA |
| ran_pars | Residual | sd_Observation | 4.127 | NA | NA | NA |
This model shows a roughly 11.1% decrease in kelp per degree increase in temperature at a site.
To further support the conclusion that temperature drives these kelp forests, we conducted a differencing analysis of those sites sampled throughout the entire duration of the study. For each year sampled at these sites, we calculated the change in S. latissima and compared it to the change in temperature using a linear mixed model with site as a random intercept and percent rocky cover as a covariate. To make our results comparable to previous analysis, rather than looking at difference in abundance per se, we looked at percent change, i.e. (kelpt - kelpt-1)/kelpt-1. This method has the advantage of removing any signal of spatial differences - be they from a gradient or just site to site variability - and focuses instead on the dynamics of temporal change. This does restrict our sample size somewhat (n=18) relative to previous analyses, but we have fewer parameters to estimate.
| Chisq | Df | Pr(>Chisq) | |
|---|---|---|---|
| TEMP_DIFF | 5.068 | 1 | 0.024 |
| SITE_PERCENT_ROCK | 0.004 | 1 | 0.952 |
| effect | group | term | estimate | std.error | df | statistic | p.value |
|---|---|---|---|---|---|---|---|
| fixed | fixed | (Intercept) | 18.456 | 79.434 | 9 | 0.232 | 0.821 |
| fixed | fixed | TEMP_DIFF | -18.733 | 8.322 | 9 | -2.251 | 0.051 |
| fixed | fixed | SITE_PERCENT_ROCK | 0.056 | 0.942 | 9 | 0.060 | 0.954 |
| ran_pars | SITE | sd_(Intercept) | 0.000 | NA | NA | NA | NA |
| ran_pars | Residual | sd_Observation | 40.961 | NA | NA | NA | NA |
Here again we see that for one degree change we have a -18.733% change in kelp.
Looking at naive simple relationships, while there is a negative correlation between winter minimum sea surface temperature, there does not appear to be a relationship between kelp and total species richness, as is found in other kelp systems.
This lack of a kelp-richness relationship could be real, or it could be that it is obscured by effects of temperature on species richness - either due to spatial or temporal variation. In other kelp systems, environmental drivers often have direct effects on bioloical communities as well as indirect effects by altering kelp abundances which then affect community structure. To tease apart the effects of temperature versus kelp on species richness, we fit a Structural Equation Model to the data using piecewise approaches with piecewiseSEM 2.0 (Lefcheck, Byrnes, and Grace 2018). In the first model, we used the same group mean centered approach for both richness and kelp abundance as above. The model was saturated, so we could not estimate fit, but, could look at
| Response | Predictor | Estimate | Std.Error | DF | Crit.Value | P.Value | |
|---|---|---|---|---|---|---|---|
| S_LATISSIMA_SQ_M + 1 | group_temp | -0.698 | 0.128 | 120 | -5.458 | 0.000 | *** |
| S_LATISSIMA_SQ_M + 1 | temp_group_cent | -0.118 | 0.048 | 120 | -2.466 | 0.014 |
|
| S_LATISSIMA_SQ_M + 1 | PERCENT_ROCK | 0.004 | 0.003 | 120 | 1.140 | 0.254 | |
| TOTAL_RICHNESS | group_temp | 0.040 | 0.100 | 120 | 0.400 | 0.689 | |
| TOTAL_RICHNESS | temp_group_cent | -0.077 | 0.015 | 120 | -5.122 | 0.000 | *** |
| TOTAL_RICHNESS | S_LATISSIMA_SQ_M | -0.002 | 0.004 | 120 | -0.660 | 0.509 | |
| TOTAL_RICHNESS | PERCENT_ROCK | 0.002 | 0.001 | 120 | 1.961 | 0.050 |
|
The relationship between temperature and richness is direct. There is no apparent mediation by sugar kelp.
Further, if we repeate the differencing analysis, but this time as an SEM with richness, we see again that temperature change, not kelp change, drives total species richness.
| Response | Predictor | Estimate | Std.Error | DF | Crit.Value | P.Value | Std.Estimate | |
|---|---|---|---|---|---|---|---|---|
| SL_PERC_DIFF | TEMP_DIFF | -18.733 | 8.322 | 4.95 | 9 | -2.251 | -0.600 | *** |
| SL_PERC_DIFF | SITE_PERCENT_ROCK | 0.056 | 0.942 | 7.72 | 9 | 0.060 | 0.016 | |
| RICH_PERC_DIFF | TEMP_DIFF | -0.119 | 0.035 | 5.90 | 8 | -3.385 | -0.858 | *** |
| RICH_PERC_DIFF | SL_PERC_DIFF | -0.001 | 0.001 | 7.94 | 8 | -0.569 | -0.144 | *** |
| RICH_PERC_DIFF | SITE_PERCENT_ROCK | -0.004 | 0.003 | 6.76 | 8 | -1.256 | -0.255 | *** |
#from 3 to 2 - line 3 #### Discussion
I. Whoah, New England is different.